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AN INVESTIGATION OF THE EFFECT OF A TURBULENT BOUND- 
ARY LAYER UPON THE LIFT OF A WING-BODY COMBINATION 
WITH A GAP BETWEEN WING AND BODY AT MACH NUMBER 1.9 

SUMMARY 

There have been several theoretical investigations which 
predict the effects of a gap between the wing and body upon 
the lift of the wing-body combination. In this investigation, 
the theoretical predictions of the lift variation with gap, 
for a given model, are corrected for the displacement thick- 
ness due to a turbulent boundary layer in the gap. These 
results are then compared with the experimental lift varia- 
tion of the model. 

The loss in lift due to a gap as determined experimentally, 
for gap widths greater than the boundary layer thickness, is 
much less than that predicted by theory. The corrections to 
the theoretical predictions made by reducing the gap by the 
displacement thickness are so small that they fail In large 
measure to explain the great difference between experiment 
and theory. 

The results of lift variation with gap for gap widths 
less than the boundary layer thickness are uncertain and this 
area of the gap problem is recommended for further investi- 
gation. 






C TO 



2 



SYMBOLS 



a 
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g 

H 

k 



L 

L g=o 

M 

n 



P 

Pt 

r 0 

3 S 

T 




V 

x 

y 

oL 

r 

S 

S 



speed of so’ind 

specific heat at constant pressure 
gap width between wing and body 
<f*/9 , boundary layer shape parameter 
thermal conductivity 

ratio of lift with gap to the lift without gap 
Mach number 

denominator of exponent describing velocity profile, 

. . . . .1/n } 

ie. u/ui- ( y/tT ) 

static pressure 

total pressure 

radius of model body 

wing body semi-span v/ith no gap 

absolute temperature 

gap parameter 

velocity component in boundary layer parallel to the 
free stream flow direction 
speed of flow 
streamwise direction 

Height above surface perpindicular to the surface 

angle of attack 

ratio of specific heats, Cp/C v 

boundary layer thickness 

boundary layer displacement thickness 



e 



boundary layer rorentmr thickness 
mass density 
absolute viscosity 



f 

r 

Subscripts 

1 refers to conditions ahead of the normal sh^ck wave 
or to the conditions at the edre of the boundary 
layer 

2 refers to conditions behind the normal shock wave 

c refers to values corrected for 

t refers to total or stagnation conditions 

oo refers to remote conditions 

1/7 

1/7 refers to the velocity profile, u/u^ = (y/aO ' 

lA.33 refers to the velocity profile, u/u^ = 

cal refers to values calculated from theory 

obs refers to values calculated from experiment 



If 



AN INVESTIGATION OF THE EFFECT OF A T TT RBULCNT BOHID- 
ARY LAYER UPON THE LIFT OF A WING- BODY CONBINATION 
••/ITH A GAP BETWEEN WING AND BODY AT I AC II NUKBBR 1.9 

INTRODUCTION 

V/lth the more frequent use of all movable control sur- 
faces on guided missies and high-performance aircraft, it has 
become of some interest to determine the effects upon lift of 
a gap between the wirg and body. There have been several 
theoretical investigations s ich as Ref. 1, to determine these 
effects. Such theoretical attempts have, of necessity, ig- 
nored the effects of a boundary layer in the gap. The present 
investigation is an attempt to determine experimentally the 

! effect of a turbulent boundary layer on a wing-body combina- 

tion with gap in supersonic flow and compare it with the the- 
oretical predictions of Ref. 1. 

The method of analysis is to first determine the bound- 

' 

ary layer displacement thickness from the measured velocity 

' 

profile. The theoretical calculation of the lift on the body 
with the measured gap versus the theoretical lift on the body 
with the gap reduced by the displacement thickness can then 
be compared with the changes in lift actually measured. 

The experimental measurements of the changes in lift 
were taken from Ref. 2, which was an investigation carried 
out simultaneously with this one at the same Nach number 



and ”ith the sa r e model in several c^nfigurati ns, including 
the crr.f i juration used in this report. 

The determination of the boundary layer displacement 
thickness in an asymmetric flow is complex and beyond the 
scope of this investigation. Of necessity, the body and wing 
were restricted to a zero degrees angle of attack. 

This investigation was conducted by LT B. J. Cartwright, 
USN, and LT Robert C. Wood, USN, during January and February, 
1955 * at the Supersonic V/ind Tunnel, University of Michigan, 
Ann Arbor, Michigan, as part of the third year curriculum in 
Aeronautical Engineering of the U. S. Naval Postgraduate 
School, Monterey, California. The project was financed by 
the Bureau of Aeronautics, Navy Department, Washington, D.C. 

The authors are especially indebted to Arnold M. Kuethe, 
Felix Pavlovski Professor of Aerodynamics, University of 
Michigan, and to H. P. Liepman, Director, Supersonic V/ind 
Tunnel, University of Michigan, for their guidance and advice 
in the course of the investigation. Grateful acknowledgement 
is also made to LT J. F. Ahearn, USN, LT K. B. Mattson, USN, 
and to the members of the V/ind Tunnel Staf T " for their help 



and assistance 
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E'WWENT AMD PR 0T-, D 1 IDE 

The experimental work was performed in t v ’e University 
of Michigan 3 inch by 13 inch Supersonic Wind Tunnel, Mach 1.9 
channel. A complete description of this facility is contained 
in Ref. 3. 

Description of Model 

The model consists of a cone cylinder with swept wings 
as shown in Fig. la. A photograph of the model mounted in 
the test section of the wind tunnel is shown in Fig. 2. 

The body has a diameter of 1.0 55 inches and an over-all 

length of 9.117 inches. The nose of the body is a cone 2.992 
inches in length with a semi-vortex angle of 10 degrees. 

The delta shaped wings measure 1.5^8 inches from tip to 
mot and have a leading edge sweep-back an^le of 35*75 
^e°:rees. To permit a gap between wing a^d body, the wings 

are attached to the body by ^eans of a pin which fits into a 

boshing in the body. In order to withstand any transient 
loads on the wing, which might occur during starting or stop- 
ping the t^rnnel, a maximum pin diameter of 0.135 inches was 
used. The arrangement is such that the leading edge of the 
wing at the root is 1.350 inches behind the shoulder of the 
nose cone. The gap between wing and body was set at O.lDD 
inches. 

To insure a turbulent boundary layer, a small wire ring 
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is f^xed to the o^se of the cone to act as a bo mdary layer 
triprino device. The ring is made of wire having a .0055 
inch diareter and is located 1.323 inches ahead of the should- 
er of the core. 

Boundary Layer Survey Equipment 

The boundary layer velocity profile was determined by 
pressure measurements with the total head probe shown in 
Fig. 1 b. The probe was constructed of stainless steel hypo- 
dermic tubing with an outside diameter of 0.042 inches. This 
tubing vas flattened and stoned to an elliptic shape with 
0.0127 by 0.0680 inches outside dimensions and 0.0063 hy 
0.0600 inches inside dimensions. With the probe against the 
body, the center of the orifice was 0.1375 mm. (0.0054 inches) 
above the body. The tubing was then soldered into 0.125 inch 
copper tubing and ground and buffed to a streamlined shape. 

The response time of the probe is six seconds. 

The boundary layer survey was made in the gap at a 
point even with the leading edge of the wing at its root. 

To traverse the boundary layer, the after part of the probe 
was securely taped to the after-body of the model with shims 
inserted beti^een the probe and the model to give the desired 
height of the probe orifice above the body. This allowed 
only one height to be used during each rim. 

This height was determined by arc measurement with a 
’Wild T-2 universal theodolite. The theodolite can normally 
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be read accurately to one second of arc. It was placed 11.4 
feet from the model, this distance being measured to a mini- 
mum accuracy of 1/8 Inches. The vibration of the probe with 
the tunnel running somewhat lessened the normal accuracy of 
the theodolite readings. Another source of error was the 
refraction of the light waves due to the density gradients 
in the boundary layer. The magnitude of this error was small, 
as indicated by measurements of the model body with and with- 
out the tunnel running. 

With these factors, the consistant accuracy of the arc 
measurements was within 8 seconds of arc, corresponding to 
0.005 inches in the probe height. 

The pressures were measured with mercury manometers which 
were read to the nearest 0.01 inches of mercury. 

Several Schlieren photographs were taken, one of which 
is shown in Fig. 3» Upon enlargement, the Schlierens showed 
the wavy fringe that typifies the edge of a turbulent bound- 
ary layer. 

Experimental procedure and Data Reduction 

Prior to each run, theodolite readings were taken of the 
body and the probe. Readings of the theodolite and manometers 
were taken during the run after the flow field and manometers 
had stabilized. 

In accordance with standard procedure for this wind 

tunnel, the remote stagnation pressure, p. 

*'00 



was taken as 
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the barometric pressure* The manometer readings ’--ere corrected 
for the difference in manometer and barometric temperatures. 

As a check on the free stream Mach number for each run, 
a reading was made of the static pressure on the floor of the 
tunnel 0.65 inches forward of the nose of the model* 

The determination of the velocity from the total head 
probe readings was made using the tables of Ref. 4. The re- 
lations used are as follows* 

For a blunt total head probe in a supersonic flow, the 
pressure measured is that behind a normal shock wave. The 
ratio of total pressures ahead of and behind the normal shock 
can be expressed as a function of the Mach number ahead of 
the shock by the following equation (see Symbols): 
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The Rayleigh pitot formula relates the ratio of the up- 
stream static pressure and the stagnation pressure behind 
the normal shock wave to Mach number as follows: 
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The ratio of velocity to the stagnation speed of sound 
is given by: 



Eq.(3) 
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These three equations are valid for the adiabatic flow 
of a perfect fluid. 

The customary additional assumptions in determining the 
velocity profile and the displacement thickness are: 

(1) The static pressure is constant through the boundary 
layer . 

(2) The Prandtl number is unity. 

(3) The thermal conductivity, k, and the absolute viscosity, 
are constant through the boundary layer. 

( 4 ) No heat is transferred through the surface of the model, 
ie., an insulated surface. 

(53 The velocity profile of the turbulent boundary layer 
obeys a power law, u/u^= (y/</")^' n . 

To these was added an additional assumption that 
(6) The loss in total pressure across the conical shock wave 
for a nose cone of a semi-vertex angle of 10 degrees is neg- 
ligible. 

By the latter assumption p]_^ r Pt^ 

This assumption was originally based on Ref. 5> which 
has condensed some of the conical shock wave data from Ref. 

6 into convenient charts and curves, one of which plotted 
the stagnation pressure ratio across conical shock waves 
versus Mach number for various cone angles. This curve, 
readable to three decimal places, indicated a negligible 
pressure loss across the shock wave for this nose cone. 

To check the validity of the assumption of a negligible 
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pressure loss across the shock waves and to obtain the local 
static pressure, three preliminary rims were made with the 
total head probe outside the boundary layer and with a needle 
static pressure probe located on the body with its orifice 
diametrically opposite to the total head probe. 

The measured static pressures were compared with the 
static pressures calculated using the barometric pressure as 



Pt l* 

Since the quantities /p^ and p^/p^. are tabulated 

2 1 2 

in Ref. 4 from equations (1) and (2), the calculation con- 
sisted of merely entering the tables with p+. /p*. and obtain- 

z 2 Z 1 

ing the corresponding p^/p^. • Then 
p i = (p i /p t 2 ) P t 2 

The calculated static pressure and the values measured 
with the static pressure probe are recorded in Table I. 

They compare closely and support the assumption of a negli- 
gible pressure loss across the shock waves. 

Since the barometer, and therefore Pt^j was seldom cons- 
tant during the investigation, the first, and if possible 
the last of the runs for each day were made with the total 
head probe in the free stream outside the boundary layer. 

The static pressure was then calculated as illustrated above. 

Since the velocity ratio, 7/a.j. given by equation (3) is 
also tabulated versus Mach number, the velocity ratio at the 
edge of the boundary layer u-ji/a^. was simultaneously determined. 
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The velocity ratio, ii/a^., in the boundary layer was ob- 
tained by entering the tables with p^/p^, being that 
from the free stream run or an average of the free stream 
runs during the day, and obtaining u/a^. directly. Then 



u 






“’T ‘ <VV 

These values are tabulated in Table I* 

Preliminary rectangular and logarithmic plots were made 
of the velocity profile to determine the boundary layer thick- 
ness, /. The edge of the boundary layer was taken where 
u/u^ m 0*994 from which 

/ m 1.2 mm. = 0.0472 inches 

The non-dimensional logarithmic plot and velocity pro- 
file are shown in Figs. 4 and 5* 

The two-dimensional displacement thickness is customar- 
ily defined from the equation 

-cT /-/ 



fr * * >1*1 ay 

or ° '' 

Bq.(4) f* ry (1 ) dy 



For steady two-dimensional flow, it is shown in Ref. 
that the energy equation reduces to 

♦ * 4?> - u a = b (k # + r ( 3y )2 

For assumptions (2), (3) and (4), a solution to the 
energy equation is 

c p T + i u 2 = c p T x + i ^ 



This may be written 
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By assumptions (1) and (5) 

Eq.(6) jo,//* : T/T 

Eq.(7) uA^ r ( y/«0 1//n 

Substituting equations (5), (6), and (7) into equation 



(4) yields 
Eq.(8) Z (f 
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The integrand of this equation was calculated for 21 
points over the interval and the integration was carried 
out graphically. 
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DISCUSSION AND RESULTS 

Boundary Layer Velocity Profile 

The equation of the velocity profile was determined to 
be 

u/i^ = (y/<0 1/n * (y/S) 1/4 - 83 

The boundary layer displacement thickness was found to 
be 

</* a 0.365 mm. s 0.0144 inches 

An adverse pressure gradient exists on the model aft of 
the shoulder due to the expansion. This pressure distribu- 
tion is shown for a similar model in Fig. 6. (Ref. 8) 

No quantitative data appeared available correlating 
boundary layer profiles with adverse pressure gradients for 
compressible flows. Perhaps this is due to the difficulty 
in accurately making such measurements without shock wave 
interference. However, R. E. Wilson, in his work on the 
growth of the turbulent boundary layer, (Refs. 9 and 10) 
for zero and favorable pressure gradients in compressible 
flow, has correlated the boundary layer shape parameter, 

H* «/*/£, with the Mach number outside the boundary layer. 

For his tabulations, Wilson chose a 1/7 power law for the 
velocity profile. This selection, he explains, was because 
a 1/7 power law has been used at subsonic speeds and has been 
found to apply for the boundary layer on a flat plate at super- 
sonic speeds for a moderate range of Reynold's numbers. 



To determine the shape parameter, the momentum thick- 
ness was calculated using the assumptions previously men- 
tioned* The equation is 

Eq. (9) e = <T ( l - gj- ) a(y/<n 
From which ° 

®l/b 83 = 0*1086 mm* = 0.00428 inches 

A 1/7 power profile is also shown in Fig* 5* Calculation 
of the displacement and momentum thicknesses for the 1/7 
profile by the methods above gave 

S* - 0.289 mm = 0*0114 inches 

© 2/7 = 0.09096 mm z 0.00358 inches 

The Mach number at the edge of the boundary layer at the 
survey station was 2.032. For this value and a 1/7 power law, 
Ref. 10 tabulates 

H1/7 = 3 * 131 

From the above calculations the shape parameters were 
h 1/7 = 3.171 

H lA.83 = 3 * 361 

This comparison of shape parameters is considered quite 
favorable. The numerical values compare closely with the 
experimental profile having the largest value. Since, other 
things being equal, the larger values of shape parameter cor- 
respond to boundary layers more nearly to the separation 
point, the effect of an adverse press’ire gradient should be 
reflected in a higher value of H. 

Also, the character of the change in the exponent, n, is 
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in a direction consistant with Clauser’s work in incompress- 
ible flow with an adverse pressure gradient. (Ref. 11) 

The probe effects and the interference effects due to 
the presence of the pin, the tape, and the reflected shock 
waves in the gap are unknown and were neglected. 

Effects of Displacement Thickness 

The authors of Ref. 1, Dugan and Hikido, in their linear- 
ized inviscid theory, have developed the variation of lift 
due to gap width as a function of a gap parameter, £ 



e + SJ ’ 

where g is the width of the gap and S* is the wing-body 
semi-span with no gap. 

The lift versus gap parameter curves were presented for 
values of r 0 /So of 0.500 and 0.216. The value of r 0 /S£ for 
the model used in this investigation is 0.25k. A curve for 
this value was derived from linear interpolation of the other 
two curves and is shown in Fig. 7. 

The corresponding experimental curve of the lift varia- 
tion with gap parameter was developed in Ref. 2. This is al- 
so shown in Fig. 7* 

Since the effect of the displacement thickness on the 
lift variation decreases as the gap increases, the curves 
are shown only for the range of gap parameter from 0 to 0.14. 
The effects of will first be considered for gap parameters 
corresponding to gap widths greater than the boundary layer 
thickness. 

The determination of a boundary layer displacement thick- 
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ness in this investigation was, of necessity, limited to two- 
dimensional flow at a zero angle of attack. The effect of 
an angle of attack would be to add a cross-flow component and 
would of course involve a three-dimensional displacement 
thickness. To evaluate these effects, however, reference was 
made to Moore 's work on compressible laminar boundary layer 
flow in Refs. 12, 13, and 14. Although Moore doesn't deal 
explicitly with a cone cylinder, his determinations of the 
displacement effect on yawed circular cones and yawed infinite 
cylinders at the plane 90 degrees from the stagnation line at 
small angles of attack would indicate a tendency for a small 
decrease in the displacement thickness. Projecting this to 
the turbulent boundary layer in this case, it was assumed 
that there would be a very small decrease, negligible for 
this report, in the displacement thickness from the two- 
dimensional case at the position of the wing for angles of 
attack up to 10 degrees, the range considered in Ref. 2. 
However, there is no experimental data to substantiate this 
assumption. 

To determine the effects of the displacement thickness 
on lift from the theoretical predictions, the value of the 
lift ratio, L/Lg.g, at a given gap parameter, was compared 
with the lift ratio for a gap reduced by the displacement 
thickness, ie., for a gap parameter of — 2 — • These 

(g -«r*) ♦ s* 

values were then compared with the experimental values, Fig. 7* 
If the only deviations of experiment from theory were due to 
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the boundary layer displacement thickness, these values would, 
of course, coincide. 

The values are tabulated in Table II. To better Illus- 
trate the comparison, Fig. 8 shows the percentage Increase 
of the lift ratio predicted by Ref. 1 due to a reduction in 
the gap width by A second curve shows the percentage 

increase of the experimental lift ratios over the theoretical 
values. 

Table II and Fig. 8 show that the displacement thickness 
has a relatively small effect on the lift variation for gap 
widths greater than the boundary layer thickness. The change 
in lift is shown to vary from 0 $ at a gap parameter of .14 
to 4.0$ at a gap parameter corresponding to the boundary layer 
thickness. The percentage of the experimental values over 
the theoretical values is shown to vary from 70.3$ to 62.0$ 
over the same range. Therefore, correcting the theoretical 
predictions for a displacement thickness accounts for only a 
small part of the difference in experiment and theory. Thus 
the concept of a boundary layer displacement thickness (or 
displacement surface In this case) fails in large degree to 
explain the big difference - in the experimental values arid 
the theoretical values of the lift variation with gap. Large 
errors in the displacement thickness calculations would still 
fail to account for the major part of this difference. For 
instance, tripling the displacement thickness at a gap para- 
meter of 0.0? would only correct the theoretical predictions 
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by 7*5% compared with a 6 7% difference between the measured 
lift ratio and that predicted by theory. The greater part 
of the difference must therefore be attributed to other fac- 
tors. Some possibilities for these other factors are the 
linearization of the equations in Ref. 1, the deviation of 
the model from the slender body requirement, other viscous 
effects not considered, such as the boundary layer on the 
surface of the wing root, the presence of the pin in the gap, 
and compressibility effects. 

Consideration was next given to the range of gap widths 
less than the boundary layer thickness. Bleviss and Struble 
in Ref. 15, and others, have indicated a belief that the 
effects of viscosity are such as to cause an inflection in 
the lift ratio versus gap curve at very small gaps so that 
the curve approaches unity somewhat asymptotically, as shown 
in Fig. 7. If the displacement surface was a real surface 
instead of an artificial one, the lift ratio would, of course, 
reach unity at a gap corresponding toe f*. However, it is 
anticipated that even at this gap width there would be some 
small flow through the gap. 

In an attempt to check the inflection in the lift ratio 
curve, additional runs were made at a gap width equal to <T*. 
The runs were made using newly installed automatic data re- 
duction equipment and the results were as follows: 

ot = 5°, h — z 1.112 J o< = 10°, — L— = 1.119 

e-o g=o 
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It is considered highly improbable that the lift ratio 
increases even for gap widths between zero and <f* • The 
accuracy of the above results is therefore questionable. The 
wind tunnel test schedule did not permit the necessary re- 
runs to check these results. This area of the gap effects 
problem is considered of high interest, however, and is 
recommended as warranting more extensive investigation. 



CONCLUSIONS AND RECOMMENDATIONS 



The following conslusions and recommendations can be 
drawn from the results of this investigation: 

1* The experimental data of Ref. 2 indicates the loss 
in lift due to a gap between wing and body to be much less 
than predicted by Ref. 1. The concept of a displacement sur 
face due to a turbulent boundary layer in the gap, fails in 
large measure to account for this difference for gap widths 
greater than the boundary layer thickness. 

2. The results of an attempt to determine the lift 
ratio at a gap width corresponding to the displacement thick 
ness are uncertain and the region of gap widths less than 
the boundary layer thickness is therefore recommended for 
further investigation* 
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TABLE I 



BOUNDARY LAYER SURVEY DATA 



RUN 

NO. 


y 

mm. 


p cal. 

in.Hg 


Pobs. 

In.Hg 


Pt2 

in.Hg 


u/ a a 


ui/a a 


u/u^ 


y/y 


1 


4.3949 


3*68 


3.70 


21.24 


- - 


1.4997 


1.0000 


3 .655 


2 


4. 4152 


3.64 


3.67 


21.15 


— 


1.5030 


1.0000 


3.668 


3 


4.3882 


3.68 


3.68 


21.26 


— 


1.5020 


1.0000 


3.<f53 




2.4006 


3.57 


— 


20.72 


- - 


1.5040 


1.0000 


2.000 


5 


0.1375 


3.57 


- - 


7.47 


0.9752 


1.5040 


0.6484 


0.114 


6 


0.1375 


3.57 


- - 


7.11 


0.9453 


1.5040 


0.6285 


0.114 


7 


0.4637 


3.57 


- - 


12.01 


1.2380 


1.5040 


0.8231 


0.386 


8 


1.4271 


3.57 


- - 


20.53 


1.5000 


1.5040 


0.9973 


1.189 


9 


0.8000 


3.57 
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FIG. 2 
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FIG. 3 

SCHLTEREII PHOTOGRAPH OF HO DEL 




29 



PIG. b 

PLOT FOR THIS DETERMINATION OF TITE PROFILE POWER 




30 



FIG. 5 

TELOCITY PROFILE 
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APPENDIX 



NOTES ON THE FEASIBILITY OF USING A 
RAKE FOR BOUNDARY LAYER MEASUREMENTS 

For total head measurements of thick boundary layers, 
such as those found on the walls of wind tunnels, etc., some 
investigators have successfully used a "rake”, consisting of 
several total head probes stacked vertically with a space be- 
tween them. 

The possibility of interference effects in supersonic 
flows and the magnification of probe effects themselves are 
obvious. However, it was felt that the economy of runs which 
could be obtained with a rake warranted investigation to see 
if the accuracy obtained would be acceptable for this invest- 
igation. Hence a rake was constructed consisting of six 
probes, each similar in construction to the single probe 
used in the investigation. The probes were stacked together 
and soldered as shown in Fig. A-l. 

Measurements were made with the rake in the boundary 
layer and free stream of the 4 inch by 4 inch wind tunnel at 
Mach number 1.9. These readings were then checked against 
measurements with a single total head probe. 

Part of this data is tabulated in Table A- I. This data 
was taken during a run with the rake centered in the free 
stream of the tunnel. It can be seen that the outer two 
probes read two inches of mercury higher than the center two, 



with the latter more closely checking with the single probe 
measurements. The boundary layer measurements likewise reg- 
istered errors in the order of 5 to 13/S* On this basis, the 
idea of using a rake was discarded in favor of the single 
probe. 



FIG. A-l 



RAKE FOR BOUNDARY LAYER MEASUREMENTS 




TABLE A- I 

RAKE AND SINGLE PROBE READINGS COMPARED IN A FREE 
STREAM AT MACH NUMBER 1.9 



RAKE PROBE 
NO. 


READING 
in. Hg 


SINGLE PROBE READING 
in. Hg 


1 


37A2 


35.57 


2 


35.93 


35.57 


3 


35.52 


35.57 


4 


35.46 


35.57 


5 


35.98 


35.57 


6 


37.43 


35.57 
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